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Abstract 

Starting from the early works by Weizsacker and Bethe about fusion cycles and energy 
conversion in stars, a brief survey of thermonuclear processes in stars leading to contem- 
porary research problems in this field is given. Special emphasis is put on the physics of 
stellar and, in particular, solar neutrinos which is at the frontline of current investigations. 
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1 Introduction 



In the 1920s Eddington formulated the hypothesis that fusion reactions between light 
elements are the energy source of the stars - a proposition that may be considered as the 
birth of the field of nuclear astrophysics [1]. It was accompanied by his pioneering work on 
stellar structure and radiative transfer, the relation between stellar mass and luminosity, 
and many other astrophysical topics. Atkinson and Houtermans [3] showed in more detail 
in 1929 - after Gamow [2] had proposed the tunnel effect - that thermonuclear reactions 
can indeed provide the energy source of the stars: they calculated the probability for 
a nuclear reaction in a gas with a Maxwellian velocity distribution. In particular, they 
considered the penetration probability of protons through the Coulomb barrier into light 
nuclei at stellar temperatures of 4- 10 7 i^. /,From the high penetration probabilities for the 
lightest elements they concluded that the build-up of alpha-particles by sequential fusion 
of protons could provide the energy source of stars. An improved formula was provided 
by Gamow and Teller [4]. 

Hence, hydrogen and helium (which were later - in the 1950s - identified as the main 
remnants of the big bang) form the basis for the synthesis of heavier elements in stars - 
but details of the delicate chain reactions that mediate these processes remained unknown 
until 1938. This is in spite of the fact that rather precise models of the late stages 
of stellar evolution existed or were soon developed. At that time, white dwarfs were 
generally considered to be the endpoints of stellar evolution, although Zwicky and Baade 
had speculated in 1934 that a neutron star could be the outcome of a supernova. In 1939 
Oppenheimer and Volkoff presented the first models of neutron stars as final stages of 
stellar evolution. Together with the so-called "frozen" or "collapsed" stars - which were 
re-named by Wheeler in 1967 as "black holes" - Chandrasekhar included these results 
for sufficiently massive progenitors in his book about stellar structure [5]. Eddington 
strongly rejected the proposal, but it proved to be true when the first rotating neutron 
star (pulsar) was detected in 1967 by Bell and Hewish. 

Probably the most important breakthrough regarding the recognition of fusion cycles 
occured in 1937/8 when Weizsacker [6,7] and Bethe [9] found the CNO-cycle - which 
was later named after their discoverers Bethe- Weizsacker-cycle (figure 3) - in completely 
independent works, and Bethe and Critchfield [8] first outlined the proton-proton chain 
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(figure 2). After a brief review of thermonuclear reactions in section 2, these nucleosyn- 
thesis mechanisms are reconsidered in section 3. 

After World War II, stellar nucleosynthesis was studied further by Fermi, Teller, 
Gamow, Peierls and others, but it turned out to be difficult to understand the forma- 
tion of elements heavier than lithium-7 because there are no stable nuclei with mass 
numbers 5 or 8. In 1946 Hoyle interpreted the iron-56 peak in the relative abundances of 
heavier elements vs. mass (figure 1) as being due to an equilibrium process inside stars 
at a temperature of 3 • 10 9 K. Later Salpeter showed that three helium nuclei could form 
carbon-12 in stars, but the process appeared to be extremely unlikely. To produce the 
observed abundances, Hoyle predicted an energy level at about 7 MeV excitation energy 
in carbon-12, which was indeed discovered experimentally, generating considerable excite- 
ment and progress in the world of astrophysics. Burbidge, Burbidge, Fowler and Hoyle 
then systematically worked out the nuclear reactions inside stars that are the basis of the 
observed abundances and summarized the field in 1957 [10]. 

The role of stellar neutrinos was considered by Bethe in [9]. Neutrinos had already 
been postulated by Pauli in 1930 to interpret the continuous beta-decay spectra, but could 
not be confirmed experimentally until 1952 by Cowan and Reines. Bethe argued in 1938 
that fast neutrinos emitted from beta-decay of lithium-4 (which would result from proton 
capture by helium-3) above an energy threshold of 1.9 MeV might produce neutrons in 
the outer layers of a star. However, this required the assumption of long-lived lithium-4, 
which turned out to be wrong. Bethe did not pursue stellar neutrinos further in his early 
works, and he or Weizsacker also did not explicitely consider the role of neutrinos in the 
initial p-p reaction, or in the CNO-cycle at that time. 

In a normal star, electron neutrinos that are generated in the central region usually 
leave the star without interactions that modify their energy. Hence, the neutrino energy is 
treated separately from the thermonuclear energy released by reactions, which undergoes 
a diffusive transport through the stellar material that is governed by the temperature gra- 
dient in the star. Stellar neutrinos are generated not only in nuclear burnings and electron 
capture, but also by purely leptonic processes such as pair annihilation or Bremsstrahlung. 

Neutrinos from nuclear processes in the interior of the sun should produce a flux of 
10 11 neutrinos per cm 2 second on the earth. In 1967 Davis et al. - following suggestions 
by Pontecorvo, and by Bahcall and Davis to use neutrinos ".. to see into the interior 
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of a star and thus verify directly the hypothesis of nuclear-energy generation in stars" 
- indeed succeeded to measure solar neutrinos with a detector based on 390000 liters of 
tetrachloroethylene. When electron- neutrinos travelling from the sun hit the chlorine-37 
nuclei, they occasionally produced argon-37 nuclei, which were extracted and counted by 
their radioactive decay. First results were published in 1968 [11]. However, these were 
neutrinos with higher energies produced in a side branch of the proton-proton chain. 
More than 90 per cent of the neutrinos are generated in the initial p-p reaction, and 
these were observed first by the Gallex collaboration in 1992 [12] using gallium-71 nuclei 
as target in a radiochemical detector. Together with the corresponding results of the 
Sage collaboration, this confirmed experimentally the early suggestions by Weizsacker 
and Bet he that p-p fusion is the source of solar energy. 

The measurements [11,12] showed that less than 50 per cent of the solar neutrinos 
that are expected to arrive on earth are actually detected. The subsequent controversy 
whether this is due to deficiencies in the solar models, or caused by flavor oscillations was 
resolved at the beginning of the 21st century by combined efforts of the SuperKamiokande 
[14] and SNO [15] -collaborations in favor of the particle-physics explanation: Neutrinos 
have a small, but finite mass, and hence, they can oscillate and therefore escape detection, 
causing the "solar neutrino deficit" - with the size of the discrepancy depending on energy. 
The identification of oscillating solar neutrinos [15] was actually preceded by evidence for 
oscillations of atmospheric muon-neutrinos - most likely to tau-neutrinos [13]. Origin of 
atmospheric neutrinos are interactions of cosmic rays with particles in the earth's upper 
atmosphere that produce pions and muons, which subsequently decay and emit electron- 
and muon-neutrinos (or antineutrinos). 

Oscillation experiments are, however, only sensitive to differences of squared masses. 
Hence, the actual value of the neutrino mass is still an open issue, and presently only the 
upper limit of the mass of the antielectron-neutrino can be deduced from tritium beta 
decay to be 2.2eV/c 2 [16]. ^From neutrinoless double beta decay [17], a lower limit of 
0.05 — 0.2eV/c 2 has been deduced in 2002, but this result is only valid if the neutrino is 
its own antiparticle, which is not certain. 

Solar neutrinos are considered in section 4, and a brief outline of some of the perspec- 
tives of the field is given in section 5. 
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2 Energy Evolution in Stars 



Stellar and in particular, solar energy is due to fusion of lighter nuclei to heavier ones, 
which is induced by thermal motion in the star. According to the mass formula that 
was derived by Weizsacker in 1935 [18], and by Bethe and Bacher independently in 1936 
[19], the difference in binding energies before and after the reaction - the mass defect 
- is converted to energy via Einstein's E = mc 2 [20], and is then added to the star's 
energy balance. The binding energy per nucleon rises with mass number starting steeply 
from hydrogen because the fraction of the surface nucleons decreases, then it flattens and 
reaches a maximum at iron-56, the most tightly bound nucleus; afterwards it drops slowly 
towards large masses. Although this smooth behavior of the fractional binding energy per 
nucleon is modified by pairing and shell effects, the overall shape of the curve ensures that 
energy can be released either by fission of heavy nuclei or by fusion of light nuclei, as it 
occurs in stars, thus providing our solar energy. 

In case of main-sequence stars such as the sun there are no rapid changes in the star 
that could compete with the time-scale of the nuclear reactions and hence, the energy 
evolution occurs through equilibrium nuclear burning. Most important in the solar case 
is hydrogen burning, where the transformation of four hydrogen nuclei into one helium-4 
nucleus is accompanied by a mass loss of 0.71 per cent of the initial masses, or 0.029m. 
It is converted into an energy of about 26.2 MeV, including the annihilation energy of 
the two positrons that are produced, and the energy that is carried away by two electron 
neutrinos. /,From the known luminosity of the sun, one can calculate a total mass loss 
rate of 4.25-10 9 A;g / ' s. At this rate, the hydrogen equivalent of one solar mass could sustain 
radiation for almost 10 11 years. 

The reactions between nuclei inside stars are due to the thermal motion, and are 
therefore called thermonuclear . Before stars reach an explosive final (supernova) stage, 
the energy release due to these reactions is rather slow. From the hydrostatic equilibrium 
condition in the sun one derives the central temperature as 

8G/iM 

With the gas constant R, the average number of atomic mass units per molecule /i (=0.5 
for ionized hydrogen), the gravitational constant G, the solar mass M = 1.99 • 10 30 kg 
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and the solar radius it!© = 6.96 • 10 8 m one finds the central solar temperature 

T c < 3 • itfK. (2) 

Numerical solutions by Bahcall et al. [22] yield a central temperature T c = 1.57- 10 7 K and 
a central pressure P c = 2.34 • 10 16 Pa, with a central solar density of p c = 1.53 ■ 10 5 kg/m 3 . 
For these large values of temperature, the assumption of an ideal gas is indeed justified. 
The reaction rates are strongly dependent on temperature (typically ~ T 22 for the CNO- 
cycle and ~ T 4 for the pp-chain at T c ) and therefore, massive stars have much greater 
luminosities with only slightly higher central temperatures. As was noted by Bethe already 
in 1938, Y Cygni has T = 3.2 • 10 7 -K" and a luminosity per mass unit of 0.l2W/kg, 
whereas the sun's luminosity per mass unit is only about 2 • 10~ 4 W/kg (the most recent 
best-estimate value [21] of the total solar luminosity being 3.842 • 10 26 W). 

Expressed in units of energy, however, the central solar temperature is only about 
1.35 ■ 10 3 keV. This has to be compared with the height of the Coulomb barrier 

Ecoui = 1 ^ — (3) 

with the interaction radius R and the proton numbers Zi, Z 2 of the nuclei that tend to 
fuse in order to release energy. Since E cou i(R) ~ Z\Z 2 MeV , more than a factor of 10 3 in 
thermal energy is missing in order to overcome the Coulomb barrier. 

Thermonuclear reactions in stars can therefore only occur due to the quantum- 
mechanical tunneling that was established by Gamow [2]. The tunneling probability 
is 

P = Po E- l ' 2 exp(-2G) (4) 

with the Gamow-factor 

_ fm 2-nZ l Z 2 e 2 

G -VY hEW ■ (5) 

Here m is the reduced mass and Zi, Z 2 are the respective charges of the fusing nuclei, 
and E is the energy. The factor p depends only on properties of the colliding system. 
For the pp-reaction at an average energy and at solar temperature, P is of the order of 
1CT 20 . It steeply increases with energy and decreases with the product of the charges. 
Hence, at solar temperatures only systems with small product of the charges may fuse, 
and for systems with larger Z 1 Z 2 the temperature has to be larger to provide a sizeable 
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penetration probability. As a consequence, clearly separated stages of different nuclear 
burnings occur during the evolution of a star in time. 

Once the Coulomb barrier has been penetrated, an excited compound nucleus is 
formed, which can afterwards decay with different probabilities into the channels that 
are allowed from the conservation laws. The energy of outgoing particles and gamma- 
rays is shared with the surroundings except for neutrinos, which leave the star without 
interactions. 

Energy levels of the decaying compound nucleus above or below the nucleon removal 
energy can be of different types, stationary levels of small width which decay via gamma- 
emission, and short-lived quasi- stationary levels above the removal energy which can 
also (and more rapidly) decay via particle emission. Their width becomes larger with 
increasing energy and eventually also larger than the distance between neighbouring levels. 

Due to the existence of quasi- stationary levels above the nucleon removal energy, a 
compound nucleus may also be formed in a resonance when the initial energy matches 
the one of an energy level in the compound nucleus. At a resonance, the cross-section can 
become very large, sometimes close to the geometrical value. Astrophysical resonant or 
non-resonant cross-sections are usually written as 



with the astrophysical cross-section factor S that contains the properties of the corre- 
sponding reaction. Although it can be computed in principle, laboratory measurements 
are a better option. However, because of the small cross-sections, these measurements 
are difficult at low energies. Extrapolations to these energies are fairly reliable for non- 
resonant reactions where S(E) is a slowly varying function of E, but this is not true in 
the case of resonances, which may (or may not) be hidden in the region of extrapolation. 
The present state of the art for measurements of S(E) in an underground laboratory to 
shield cosmic rays is shown in figure 4 for the reaction 



that is very important in the stellar pp-chain, cf. next section. The solid line is a fit with 
a screening potential that accounts for a partial shielding of the Coulomb potential of 
the nuclei due to neighbouring electrons. Data from the LUNA collaboration [23] extend 



a(E) = S ■ E^expi-ZG) 



(6) 




(7) 
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down to 21 keV, where the Gamow peak at the solar central temperature is shown in 
arbitrary units. The peak arises from the product of the Maxwell distribution at a given 
temperature T and the penetration probability. Its maximum is at an energy 

r fm 27rZ 1 Z 2 e 2 A;T 1 2/3 
Eg=[^ ~ h j • (8) 

At Eq, the S-factor for the He-3 + He-3 reaction becomes 5.3MeVb. The average 
reaction probability per pair and second is given by 

<av>= / a(E)vf{E)dE (9) 
Jo 

where f(E) can be expressed in a series expansion near the maximum. Keeping only the 
quadratic terms, the reaction probability becomes [24] 

<™>=i&"w) ! * Se -^ a *- T) (10) 

with the S-factor Sq at the Gamow peak and 

r = 3E G /(kT). (11) 



The temperature dependence of < av > may be expressed as 

din < av > t 2 
dtnT ~ 3 ~ 3' 



;i2) 



which can attain values near or above 20. As a consequence of such large values for the 
exponent of T, the thermonuclear reaction rates become extremely strongly dependent on 
temperature, and small fluctuations in T may cause dramatic changes in the energy (and 
neutrino) production of a star. The corresponding uncertainty in stellar models created 
the long-standing controversy about the origin of the solar neutrino deficit, which has 
only recently been decided in favor of the particle-physics explanation, cf. section 4. 
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3 Hydrogen burning 



Due to the properties of the thermonuclear reaction rates, different fusion reactions in 
a star are separated by sizeable temperature differences and during a certain phase of 
stellar evolution, only few reactions occur with appreciable rates. Stellar models account 
in network-calculations for all simultaneously occuring reactions. Often the rate of the 
fusion process is determined by the slowest in a chain of subsequent reactions, such as in 
case of the nitrogen-14 reaction of the CNO-cycle. 

In hydrogen burning, four hydrogen nuclei are fused into one helium-4 nucleus, and the 
mass defect of 0.71 per cent is converted into energy (including the annihilation energy 
of the two positrons, and the energy carried away by the neutrinos): 

4- 1 H ^ 4 He + 2e + + 2u e + 26.2MeV. (13) 

As net result, two protons are converted into neutrons through positron emission {beta + - 
decay), and because of lepton number conservation, two electron neutrinos are emitted. 
Depending on the reaction which produces the neutrinos, they can carry between 2 and 
30 per cent of the energy. Helium synthesis in stars proceeds through different reaction 
chains which occur simultaneously. The main series of reactions are the proton-proton 
chain, figure 2, and the CNO-cycle, figure 3. 

In the present epoch, the pp-chain turns out to be most important for the sun - the 
CNO-cycle produces only 1.5 per cent of the luminosity [22]. The pp-chain starts with 
two protons that form a deuterium nucleus, releasing a positron and an electron neutrino. 
(With much smaller probability it may also start with the p-e-p process, figure 2). This 
reaction has a very small cross section, because the beta-decay is governed by the weak 
interaction. At central solar temperature and density, the mean reaction time is 10 10 
years, and to a certain extent it is due to this huge time constant that the sun is still 
shining. With another proton, deuterium then reacts to form helium-3. This process is 
comparably fast and hence, the abundance of deuterons in stars is low. 

To complete the chain to helium-4 three branches are possible. The first - in the sun 
with 85 per cent most frequent - chain (ppl) requires two helium-3 nuclei and hence, the 
first reaction has to occur twice, with two positrons and two electron neutrinos being 
emitted. The other two branches (ppll, ppIII) need helium-4 to be produced already 
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(in previous burnings, or primordially) . In the subsequent reactions between helium-3 
and helium-4, the additional branching occurs because the product beryllium-7 can react 
either with an electron to form lithium-7 plus neutrino (ppll), or with hydrogen to form 
boron-8 (ppIII). The energy released by the three chains differs because the neutrinos 
carry different amounts of energy with them, and the relative frequency of the different 
branches depend on temperature, density, and chemical composition. The per centages in 
figure 2 refer to the standard solar model at the present epoch [22]. Details of the various 
parts of the chain including the corresponding energy release, the energies carried away 
by the neutrinos and the reaction rate constants have been discussed by Parker et al. [26] 
and Fowler et al. [27]. 

The other main reaction chain in hydrogen burning is the CNO-cycle, figure 3. Here, 
the carbon, nitrogen and oxygen isotopes serve as catalysts, their presence is required for 
the cycle to proceed. The main cycle is completed once the initial carbon-12 is reproduced 
by nitrogen- 15 + hydrogen. There is also a secondary cycle (not shown in figure 3 since 
it is 10 4 times less probable). It causes oxygen-16 nuclei which are present in the stellar 
matter to take part in the CNO-cycle through a transformation into nitrogen-14. The 
CNO-cycle produces probably most of the nitrogen-14 found in nature. For sufficiently 
high temperatures, the nuclei attain their equilibrium abundances and hence, the slowest 
reaction - which is nitrogen-14 + hydrogen - determines the time to complete the whole 
circle (bottom of figure 3). 

The CNO-cycle contributes only a few per cent to the luminosity of a star with one 
solar mass, but it dominates in stars with masses above 1.5 times the solar value because 
its reaction rates rise much faster with temperature as compared to pp. Details of the 
Bethe-Weizsacker cycle have been discussed by Caughlan and Fowler [28]. The cycle had 
first been proposed by Weizsacker in [7]. In this work, he abandoned the main reaction 
path that he had considered in [6], namely, from hydrogen via deuterium and lithium to 
helium, because the intermediate nuclei of mass number 5 that were supposed to be part 
of the scheme had turned out to be unstable. 

In the first paper of the series [6], he had considered various reaction chains that allow 
for a continuous generation of energy from the mass defect, and also of neutrons for the 
buildup of heavy elements. He had confirmed that the temperatures in the interior of stars 
are sufficient to induce nuclear reactions starting from hydrogen. In the second paper he 
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modified the results; in particular, he discussed the possibility that some of the elements 
might have been produced before star formation by another process. 

The link between energy evolution in stars and the formation of heavy elements as con- 
sidered in [6] turned out to end up in difficulties when calculated quantitatively. Hence, 
he modified his version of the so-called " Aufbauhypothese" , according to which the neu- 
trons necessary for the production of heavy elements should be generated together with 
the energy, and decoupled the generation of energy from the production of heavy ele- 
ments. He then concluded that stellar energy production should essentially be due to 
reactions between light nuclei, with the corresponding abundances being in agreement 
with observations. The CNO-cycle was considered to be the most probable path. 

In his independent and parallel development of the CNO-cycle that was published 
somewhat later [9] and contained detailed calculations, Bethe showed that "... there will 
be no appreciable change in the abundance of elements heavier than helium during the 
evolution of the star but only a transmutation of hydrogen into helium. This result... is in 
contrast to the commonly accepted 'Aufbauhypothese'" . Here, he referred to Weizsacker's 
first hypothesis [6] which had, however, already been modified [7]. 

Together with Critchfield [8], Bethe also investigated essential parts of the pp-chain 
(which Weizsacker also mentioned) - in particular, deuteron formation by proton combi- 
nation as the first step - and came to the conclusion that it "...gives an energy evolution 
of the right order of magnitude for the sun" . Details of the pp-chain were developed much 
later in the 1950s by Salpeter [29] and others. In 1938/9, however, Bethe was convinced 
that "... the reaction between two protons, while possible, is rather slow and will therefore 
be much less important in ordinary stars than the cycle (1)" namely, the CNO-cycle. 

In a calculation of the energy production by pp-chain versus CNO-cycle (figure 5), 
Bethe obtained qualitatively the preponderance of H+H at low and N+H at high temper- 
atures. However, the result had to be modified in the course of time as it became evident 
that the pp-chain is more important than the CNO-cycle at solar conditions, although 
the Bethe- Weizsacker-fraction will increase considerably in the coming 4 billion years, and 
eventually supersede the contribution from the ppll-chain (figure 6). 

Today, detailed solar models allow to calculate the fractions of the solar luminosity 
that are produced by different nuclear fusion reactions very precisely [22]. The model 
results not only agree with one another - in the neutrino flux predictions to within about 
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1 per cent - they are also consistent with precise p-mode helioseismological observations 
of the sun's outer radiative zone and convective zone [30]. Moreover, the production of 
heavier elements up to iron in subsequent burnings at higher temperatures [27], as well 
as beyond iron in the r- and s-process is rather well-understood [31]. 

4 Stellar Neutrinos 

In stellar interiors, only electron neutrinos play a role. The interaction of neutrinos with 
matter is extremely small, with a cross-section of 

a v ~ (E v /m e c 2 ) 2 ■ l(T 17 m&. (14) 

Hence, the cross-section for neutrinos with E u ~ lMeV is o v ~ 3.8 • 10~ 17 mb, which is 
smaller than the cross-section for the electromagnetic interaction between photons and 
matter by a factor of about 1CT 18 . Associated with the cross-section is a mean free path 

\ v = ~ m (15) 

p-o v p 

with the atomic mass unit u = 1.66 • 10~ 27 kg and p in kg/m 3 . In stellar matter with 
p ~ 1.5 • 10 3 kg/m 3 , the mean free path of neutrinos is therefore approximately 

A„ ~ 3 • 10 17 m ~ lOpc ~ 4 • 10 9 R Q (16) 

and hence, neutrinos leave normal stars without interactions that modify their energy. 
This is different during the collapse and supernova explosion in the final stages of the 
evolution of a star where nuclear density can be reached, p ~ 2.7 ■ 10 17 kg/m 3 such that 
the mean free path for neutrinos is only several kilometers, and a transport equation for 
neutrino energy has to be applied. 

Here only the neutrinos from nuclear reactions in a normal main-sequence star like the 
sun are considered; their energies are (to some extent, since the continuous distributions 
overlap) characteristic for specific nuclear burnings. The pp-chain which provides most of 
the sun's thermonuclear energy produces continuum neutrinos in the reactions ([32]; cf. 
figure 2) 

1 H + 1 H H + e+ + v e (0A20MeV) 
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8 B ^ 8 Be* + e+ + v e (U.OQMeV) 

13 N C + e+ + z/ e (1.20MeTO 

15 iV + e+ + z/ e (1.74Me\/) 
where the numbers are the maximum neutrino energies for the corresponding reaction. 
In addition to these continuum neutrinos, there are neutrinos at discrete energies from 
the pp-chain 

1 H + 1 H + e~ H + v e (l.UMeV) 
7 Be + e~ -* 7 Li* + v e (0.861MeV - QOpercent) 

(0.383Me\/ - lOpercent) 
(depending on whether lithium-7 is in the ground state, or in an exited state) 

8 B + e" ^ 8 Be + v e {IbMMeV). 
The CNO-cycle (figure 3) which becomes important in stars with masses above 1.5 
solar masses, or in later stages of the stellar evolution (figure 8) also produces neutrinos 
at discrete energies 

13 N + e~ ^ 13 C + v e (2.22MeV) 
15 + e" N + v e (2.7QMeV). 
For experiments to detect these neutrinos when they arrive on earth 8.3 minutes after 
their creation the flux at the earth's surface is of interest. Neutrinos from the central 
region of the sun yield a flux of about 10 7 /(m 2 ■ s). The precise value as function of the 
neutrino energy can be calculated from solar models ([25], figure 7). Here, solid lines 
denote the pp-chain and broken lines the CNO-cycle. The low-energy neutrinos from the 
initial pp-reaction yield the largest flux. However, the first experiment by Davis et al. [11] 
that detected solar neutrinos on earth with a large-scale underground tetrachloroethylene 
tank in 1967/8 - and thus confirmed the theory how the sun shines and stars evolve - 
made use of the reaction 

u e +f 7 CI -> e" +H Ar - 0.8UMeV 
and hence, only neutrinos with energies above 0.814 MeV could be observed through 
the decay of radioactive argon nuclei - which are mostly the solar boron-8 neutrinos, 
cf. figure 7. The rate of neutrino captures is measured in solar neutrino units; 1 SNU 
corresponds to 1CT 36 captures per second and target nucleus. Experimental runs by Davis 
et al. during the 1970s, 80s and 90s yielded a signal (after subtraction of the cosmic-ray 
background 1.6 kilometers underground) of (2.3±0.3)SNU, whereas the predicted capture 
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rates from a solar model were SNU for pp (because it is below threshold), 5.9 SNU 
for boron-8 beta decay, 0.2 SNU for the pep- reaction, 1.2 SNU for beryllium-7 electron 
capture, 0.1 SNU for nitrogen-13 decay and 0.4 SNU for oxygen-15 decay, totally about 
8 SNU. 

The observation of less than 50 per cent of the expected neutrino flux created a 
controversy about the origin of the deficit, which was finally - in 2001 - resolved [15] in 
favor of the particle-physics explanation that had originally been proposed by Pontecorvo 
in 1968 [33]: on their way from the solar interior to the earth, electron-neutrinos oscillate 
to different flavors which escape detection, thus creating the deficit. Although deficiencies 
in the solar models could have been responsible for the discrepancy (in view of the sensitive 
dependence of the neutrino flux on the central temperature), it could be confirmed [15] 
that the models are essentially correct, giving the right value of T c within 1 per cent. 

Before this big step in the understanding of stellar evolution and neutrino properties 
could be taken, there was substantial progress both in experimental and theoretical neu- 
trino physics. In 1992 the Gallex-collaboration succeeded to measure the pp- neutrinos 
from the initial fusion reaction, which contributes more than 90 per cent of the inte- 
gral solar neutrino flux [12]. They used a radiochemical detector with gallium as target, 
exploiting the reaction 

v e +l\ Ga -> e" +H Ge - 0.23MeU. 

The threshold is below the maximum neutrino energy for pp-neutrinos of 0.42 MeV 
and a large fraction of the pp-neutrinos can therefore be detected in addition to the pep-, 
beryllium-7 and boron-8 neutrinos. Gallex - which is sensitive to electron neutrinos only - 
thus provided the proof that pp- fusion is indeed the main source of solar energy. The result 
[(69.7 + 7.8/ - 8.1) SNU] was confirmed by the Sage experiment [(69 ± 12) SNU] in the 
Caukasus [35]. Again this was substantially below the range that various standard solar 
models predicted (120-140 SNU), and the solar neutrino deficit persisted. At that time, 
there were clear indications - but no definite evidence yet - that the flux decreases between 
sun and earth due to neutrino flavor oscillations - most probably enhanced through the 
MSW-effect [36] in the sun -, "...pointing towards a muon-neutrino mass of about 0.003 
eV" [34]. The result was later updated to (73.9 ± Q.2)SNU and could be assigned to 
the fundamental low-energy neutrinos from the pp and pep reactions - but then there 
remained no room to accomodate the beryllium-7 and boron-8 neutrinos. 
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Solar neutrinos were also detected in real time with the Kamiokande detector [37] 
in Japan, a water- Cerenkov detector, and precursor to the famous SuperKamiokande 
detector. Due to the high threshold of about 7.5 MeV, it could see only the most energetic 
neutrinos from the decay of boron-8 in the solar center. With the Cerenkov light pattern 
one could measure for the first time the incident direction of the scattering neutrinos, and 
prove that they do indeed come from the sun. The result of the boron-8 neutrino flux was 
/^f=(^) = 0-54 ±0.07, 

again confirming the deficit. To solve the solar neutrino problems, a larger target 
volume and a lower energy threshold was needed: the SuperKamiokande detector in the 
same zinc mine with a threshold of 5 MeV. Here, 32000 tons of pure water are surrounded 
by 11200 photomultiplier tubes for observing electrons scattered by neutrinos (many of 
the tubes were destroyed end of 2001 when one collapsed, emitting an underwater shock 
wave). The detector was designed to record about 10000 solar neutrino collisions per year 
- 80 times the rate of its predecessors -, but also atmospheric neutrinos, and possible signs 
for proton decay. The result [14] for the boron-8 flux can be expressed as 

/i«*^^(i/ e ) = 0.47 ±0.02, 
which was in agreement with the previous findings, but more precise. There was a mas- 
sive hint that the deficit could be due to neutrino oscillations, since the SuperKamiokande 
collaboration found evidence in 1998 [13] that muon neutrinos which are produced in the 
upper atmosphere by pion and muon decays change their type when they travel distances 
of the order of the earth's radius due to oscillations into another species, most likely into 
tau neutrinos. The appearance of the tau neutrinos could not yet be detected directly, 
but oscillations to electron neutrinos in the given parameter range were excluded since 
the z/ e -flux was unchanged, and also by reactor data. Accelerator experiments with a long 
baseline of 700 km between neutrino source and detector are being planned to verify this 
interpretation [40]. 

The atmospheric data showed a significant suppression of the observed number of 
muon neutrinos as compared to the theoretical expectation at large values of x/E u , with 
the travel distance x (large when the neutrinos travel through the earth) and the neutrino 
energy E u which is in the GeV-range for atmospheric neutrinos and hence, much higher 
than in the solar case. The observed dependence on distance is expected from the the- 
oretical expression for oscillations into another flavor, which yields in the model case of 
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two flavors 




(17) 



Here, 9 is the mixing angle between the two flavors considered, and the characteristic 
oscillation length (the distance at which the initial flavor content appears again) is 



with the difference Am 2 =| m 2 , — m\ | of the neutrino mass eigenstates. Atmospheric 
neutrino experiments are thus sensitive to differences in the squared masses of 10 -4 to 
10~ 2 eV 2 /c A whereas solar neutrinos are sensitive to differences below 2- !0~ A eV 2 /c 4 due to 
the lower neutrino energy and the larger distance between source and detector. Whereas 
such mixings between neutral particles that carry mass had been firmly established many 
years ago in the case of quarks that build up the K° and _B°-mesons and their antiparticles 
- including the proof that CP is violated [38] for three quark families -, it remained an 
open question until 1998 whether the corresponding phenomenon [39] exists for leptons. 

The atmospheric SuperKamiokande data proved beyond reasonable doubt the exis- 
tence of oscillations and finite neutrino masses [13, 40], with Am 2 tm ~ 2.5 • !0~ 3 eV 2 and 
maximal mixing sin 2 {26 atm ) ~ 1. The corresponding step for solar neutrinos followed in 
2001: charged-current results [y e + d — * e + p + p; sensitive to electron neutrinos only) 
from the Sudbury Neutrino Observatory (SNO) in Canada with a -D 2 0-Cerenkovdetector 
[15], combined with elastic scattering data from SuperK (z/ + e — > v + e; sensitive to all 
flavors), established oscillations of solar boron-8 neutrinos. 

These results were confirmed and improved (5.3<x) in 2002 by neutral current results 
from SNO (v & +d — > u+n+p); a further improved measurement with salt (NaCl; chlorine- 
35 has a high n-capture efficiency) is currently underway. The total flux measured with 
the NC reaction is (5.09 + 0.64/ — 0.61) ■ 10 6 neutrinos per (cm 2 ■ s). This is in excellent 
agreement with the value from solar models (5.05 + 1.01/ — 0.81), proving that stellar 
structure and evolution is now well-understood. The currently most-favored mechanism 
for solar neutrino conversion to myon- and tauon-flavors is the "Large mixing angle" 
solution, which also implies matter-enhanced (resonant) mixing in the interior of the sun 
through the MSW-effect [36]. 



L = A-KE U /Am 2 



2A8(E u /MeV)/(c 4 Am 2 /eV 2 ) [m] 



(18) 
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5 Perspectives 



Since the early works by Weizsacker and Bethe, the investigation of thermonuclear pro- 
cesses in stars has developed into considerable detail, and with the advent of stellar 
neutrino physics an independent confirmation of the origins of solar energy has emerged. 
Improvements in the precise measurements of all the reaction rates at low energies that 
are involved in the fusion chains may still be expected, as has been outlined in the model 
case of the 3 He + 3 He system within the energy region of the solar Gamow peak [23]. 

However, not only a good knowledge of the processes involved in equilibrium burnings 
at energies far below the Coulomb barrier, but also of explosive burning (with short-lived 
nuclides at energies near the Coulomb barrier) is of interest, because both contribute to 
the observed abundances of the elements. This requires new experimental facilities. 

An improved understanding of the cross-sections will then put the predictions of the 
solar neutrino flux and spectrum on a better basis. This entire spectrum will be investi- 
gated with high precision in the coming decade. In particular, the new detector Borexino 
will measure the monoenergetic beryllium-7 neutrinos at 862 keV, which depend very 
sensitively on the oscillation parameters. The LENS experiment will utilize inverse beta- 
decay to an isomeric state of the daughter nuclide in order to investigate the low-energy 
solar neutrinos in real time with considerably reduced background. Together with forth- 
coming SNO and KamLAND results it will then be possible to definitely determine all 
the mixing parameters in a three-family scheme - and verify, or falsify, the LMA solution. 

The more detailed knowledge about the physics of stars will thus be supplemented 
by considerable progress regarding neutrino properties [40]. Questions to be settled are 
the individual neutrino masses (rather than the difference of their squares); whether 
neutrinos are their own antiparticles (to be decided from the existence or non-existence of 
neutrinoless double beta-decay); whether neutrinos violate CP just as quarks do, or maybe 
in a different manner that opens up a better understanding of the matter-antimatter 
asymmetry of the universe than has been possible from the investigation of quark systems 
(to be decided in experiments with strong neutrino beams). 

In any case, the Standard model of particle physics has to accomodate finite neutrino 
masses, and in future theoretical formulations the relation between quark mixing and 
neutrino mixing will probably become more transparent. 
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Figure captions 

Fig. 1. Solar system abundances of the nuclides relative to silicon (= 10 6 ) plotted as 
function of mass number. The stellar nuclear processes which produce the charac- 
teristic features are outlined. The p-p-chain and the CNO-cycle of hydrogen burning 
are discussed in the text. Elements up to A < 60 are produced in subsequent burn- 
ings at higher temperatures, beyond A = 60 in supernovae through the r-, s- and 
p-pro cesses. 

Source of the data: 

A.G.W Cameron, Space Sci. Rev. 15, 121 (1973). 
Source of the Figure: 

K.R. Lang, Astrophysical Formulae, p. 419. Springer (1980). 

Fig. 2. Proton-proton reactions are the main source of stellar energy in stars with masses 
close to or below the solar value. They were already briefly considered by Weizsacker 
[6] and discussed in more detail by Critchfield and Bethe [8]. Today it is known that 
the ppl branch is supplemented by the ppll branch, and the small, very temperature- 
dependent ppIII branch. In the latter two branches, additional electron neutrinos 
of fairly high energy are produced. The approximate partitions refer to the sun. 

Source of the Figure: 

H. Karttunen et al., Fundamental Astronomy. Springer (1987). 

Fig. 3. At temperatures above 20 Million Kelvin corresponding to stars of more than 1.5 
solar masses the Bethe- Weizsacker-cycle is more important than the proton-proton 
chain because its reaction rate rises faster with temperature. This CNO-cycle was 
first proposed by Weizsacker [7] and Bethe [9]. Here, carbon, oxygen and nitrogen 
act as catalysts. 

Source of the Figure: 

H. Karttunen et al., Fundamental Astronomy. Springer (1987). 



21 



Fig. 4. The astrophysical cross-section factor S(E) for the reaction 3 He( 3 He, 2p) A He. 
The solid line is a fit with a screening potential. Data from the LUNA collaboration 
[23] extend down to 21 keV, where the Gamow peak at the solar central temperature 
is shown in arbitrary units. 

Source of the Figure: 

E.C. Adelberger et al., Rev. Mod. Phys. 70, 1265 (1998). 

Fig. 5. Stellar energy production in 10~ 4 J/ (kg-s) due to the proton- proton chain (curve 
H+H) and the CNO-cycle (N+H), and total energy production (solid curve) caused 
by both chains. According to this calculation by Bethe in 1938 [9], the CNO- 
cycle dominates at higher than solar temperatures. Its role at and below solar 
temperatures as compared to pp is, however, overestimated, cf. fig. 6. 

Source of the Figure: 

H.A. Bethe, Phys. Rev. 55, 434 (1938). 

Fig. 6. Fractions of the solar luminosity produced by different nuclear fusion reactions 
versus solar age, with the present age marked by an arrow (Bahcall et al. 2001 
[22]). The proton-proton chain is seen to generate the largest luminosity fractions 
- in particular, through the branch that is terminated by the 3 He — 3 He reaction. 
The solid curve shows the luminosity generated by the CNO-cycle, which increases 
with time, but is only a small contribution today. 

Source of the Figure: 

J.N. Bahcall, M.H. Pinsonneault and S. Basu, Astrophys. J. 555, 990 (2001). 
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Fig. 7. Spectrum of solar electron neutrinos according to the "Standard Solar Model". 
The largest contribution is generated by low-energy neutrinos from the p-p chain 
that have been detected by the Gallex experiment [12,34]. Solid lines indicate neu- 
trinos from the pp-chain, dashed lines from the CNO-cycle. Neutrinos of higher 
energies had first been observed by Davis et al. [11]. The calculation is by 
Bahcall and Pinsonneault [25]. (The hep-neutrinos arise from the ppIV-reaction 
3 He + p — > 4 He + v e + e + which is not shown in figure 2). 

Source of the Figure: 

J. Bahcall and M.H. Pinsonneault, Rev. Mod. Phys. 67, 1 (1995) and |astro-ph 
001034a 



Fig. 8. The proton-proton, beryllium-7, boron-8 and nitrogen-13 neutrino fluxes as func- 
tions of solar age, with the present age marked by an arrow (Bahcall et al. [22]). The 
Standard Solar Model ratios of the fluxes are divided by their values at 4.57 • 10 9 y, 
the present solar age. 

Source of the Figure: 

J.N. Bahcall, M.H. Pinsonneault and S. Basu, Astrophys. J. 555, 990 (2001). 
[this figure should be omitted in case of space problems], 
((please obtain permission to reprint the figures)) 
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